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We present recent progresses in radiative penguin decays of B meson using a large sample of BB pairs recorded 
at the T{4S) resonance with the Belle detector at the KEKB asymmetric energy e+e~ collider. We report precise 
measurement of inclusive 6 ^ S7 branching ratio with cut Ej > 1.7 GeV, first measurement of time-dependent CP- 
violation in B'^ —* Ksp'-^-f, measurement of i3+ — > K^ri'^y branching fraction, and improved branching fraction results 
for (p, a;)7 with new CP and isospin violation results in the mode. 



1. Introduction 



O ' Radiative B decay proceeds through a penguin loop diagram in the Standard Model (SM). New particles beyond 
the SM may contribute in the loop diagrams. In this paper, we present recent Belle results on various radiative 
t — B decay topics below. Inclusive 6 — > 57 measurements will constrain New Physics from the branching fraction as 
' ' well as determining 6-quark mass and motion from photon energy spectrum. Time dependent CP asymmetry in 
1— I & — > 57 transition can probe for right-handed coupling which is not present in SM. Measurements of exclusive mode 
^ B rj' and B — *■ (p, 0^)7 can probe deviations from SM by the branching fraction and CP asymmetry. 



In all these studies, large continuum background from e+e qq events {q = u, d, s, c) has been suppressed using 



^ ' event topology variables. In exclusive modes, two kinematic variable Mtc = \/ ^-bcam ^ "Pb = ^ -^'bcam 

i are commonly used. Here i^beam is the beam energy and E*^ and •p*Q is the energy and momentum of a B candidate, 

\ all defined in the T(4S') frame. 
^ ■ 
> 

\Q 2. Inclusive sj measurement 

O ; 

. In this fully inclusive method, b ^ sj decay is studied using a 604.6fb~^ ON data sample taken at the T(4S') 

' resonance, and 68.3fb^^ OFF data sample taken at an energy 60MeV below the resonance. Hard photon is selected 

QQ ' from well isolated ECL clusters with the shower shape consistent with a photon, and required to have a center-of-mass 

' energy K > 1.4 GeV. tt" and rj vetoes are then applied. 

* * I 

J> , First, the photon energy spectrum measured in OFF data is scaled by luminosity to the expected number of 
non-BB events in ON data and subtracted, shown in Fig. [Ija). Here the OFF spectrum has been corrected due to 

^ ■ the slight lower energy than ON events. Then the backgrounds coming from non-primary photons from B mesons 
. 1 subtracted from the obtained spectrum. From Monte Carlo (MC) study, six background categories are considered 
(Fig. IHb)): (i) photons from 7r° 77; (ii) photons from rj — + 77; (iii) other real photons; (iv) ECL clusters not 
due to single photons (mainly A'^'s and n's); (v) Electrons misidentified as photons; and (vi) beam background. For 
each of the background category, the shape and yield is corrected from MC using data by selecting them. Then, 
for each selection criterion used in this analysis, the data-MC efficiency ratio is obtained using appropriate control 
samples, and is then used to scale the MC background sample. For example, the 7r° veto efficiency is studied using 
£)*+ Z3''(/f ^7r+7r°)7r+ decays where only one 7 from 7r° is reconstructed. 

The raw spectrum shown in Fig. [Ijc) is obtained after subtracting the six background categories after scaling. The 
raw spectrum is then corrected by three step procedure: (i) divide by the efficiency of the selection; (ii) perform an 
unfolding which removes the distortion by caused ECL; (iii) divide by efficiency of detection. Then, two additional 
corrections are appfied: (i) to remove cabibbo suppressed B Xdj decays, using the ratio of branching fractions; 
(ii) to correct the measurements to _B-meson rest frame from T(4S') frame (boost correction). The partiaf branching 
fraction, first moment (mean) and second central moment (variance) of the photon energy spectrum from B decay 
are obtained as: ^ B (B ^ Xsj) = (3.31 ± 0.19 ± 0.37 ± 0.01) x lO"''; (E^) = 2.281 ± 0.032 ± 0.053 ± 0.002 GeV; 
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Figure 1: (a) ON data (open circle), scaled OFF data (open square) and continuum background subtracted (filled circle) 
photon energy spectra of candidates in the c.m.s frame, (b) The spectra of photons from _B-meson decays passing selection 
criteria as predicted using a MC sample, (c) The extracted photon energy spectrum of i? — > Xs_d7. The two error bars show 
the statistical and total errors. 
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Figure 2: (Left) Mbc distributions for Kg-K'^ir^'y events in region. The dashed and dash-dotted curves are the qq and 

all BG. The thin curve is the total signal and the thick curve is the total PDF. (Right) Fit projections on the At distributions 
for events good-tagged as (a) and (b) . The raw asymmetry as a function of At is shown in (c). 



{E'^) ~ {E^f = 0.0396 ± 0.0156 ± 0.0214 ± 0.0012 GeV^, where the errors are statistical, systematic and from the 
boost correction, respectively. 



3. Time-dependent CP Asymmetries in S° K^p^i Decay 

In this mode, signal decay vertex can be reconstructed from two charged pions to calculate the decay time 
difference At between signal and tag side B, thus avoiding A'^ vertexing. The 5° Kgp'^j candidates are selected 
from the Kgir^ir^j sample by requiring the tt+tt^ invariant mass to lie in the region, 0.6GeV/c^ < mTTTr < 
0.9GeV/c'^. We first measure the effective CP-violating parameters, Scs and ^off, using the final sample and then 
convert them to the CP- violating parameters of P" Rgp'^-f using a dilution factor V. 

Fig-IHleft) shows a fit to A/bc in AP signal region after vertexing. We obtain 212 ± 17 total signal yield from the 
total 299 events in the signal Afbc region. The Scs and AcS parameters are then extracted from an unbinned maximum 
likelihood fit to the At distribution, as 5off = 0.09 ± 0.27(stat.)to °7(syst.) and Aff = 0.05 ± 0.18(stat.) ± 0.06(syst.). 
The At projection for the fit is shown in Fig. [^fright). 
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Figure 3: Projections from the 2D fit to data, (a) and (b) plots show M^c distribution in the signal region of AE and AE 
distribution in the signal region of Mhc, for the B"*" — > K'^rj'^ mode. The Kri''y function is shown in dashed red, qq in dotted 
black, 6 ^ c in dash-dotted green, 6 — > m, d, s in solid magenta, and the combined function in solid blue, (c) and (d) are the 
same for Kgrj'j mode. 



Mode 


Yield(events) 


e n 


BF(10"'^) 


S'(a) UL(1Q-^) 


B+ 




QO R + ii s 
c, 1+5.0 


0.027 0.571 
0.016 0.197 


q O+1.2+0.3 
■^■^-1.1-0.3 
r, .+2.4+0.4 
^■^-1.9-0.5 


3.3 

1.3 6.3 



Table I: The yields, efiiciencies(£), daughter BFs (n)i measured BFs, fit significances including systematics (S') and upper 
limits (UL) for the measured decays. 

The parameter Scs is related to S for Kgp'~'j with a dilution factor V: 

5eff _ f[\FA\' + 2U{FXFB) + F*s{K)FB{K)] 
5^0,0, J[\FA\^ + 2^iFXFB) + \FB\^] ' 

where Fa, Fb are photon-helicity averaged amplitudes for — > iir^p°(7r+7r^)7 and K*^ {Kgn^)Tr^ j , respec- 

tively. From the study to (i4r+7r~7r+) system in charged mode — + K^tt^tt^j using isospin symmetry, we obtain 
V = O.SSt'^ Q^. In summary, we measure the CP-violating parameter Sj^OpO^ = 0.11 ±0.33(stat.)lQ gg(syst.) [2]. This 
agrees with SM prediction Si^OpO^ « 0.03 

4. Evidence for B r]'K^ Decays at Belle 

The exclusive mode B rj' Kj is analyzed using 604.6fb~^ of data, for both charged and neutral mode. ?/ mesons 
are reconstructed as 77' t^tt+tt^ and 77' p'^j. 77 mesons are reconstructed as 77 —> 77 and 77 Tr+Tr^Tr*^. D'^ mass 
veto on mx+TT- for any changed pion in the event and J/ip mass veto on to,,'-,, are applied to suppress b —>■ c events. 

The signal yield is extracted by fitting 2D Mbc and AE distributions shown in Fig.[31 All final states are combined 
in to one distribution for fitting for the charged mode, so does the neutral mode. The signal shape is calibrated 
using large samples oi B ^ if* (892)7 data and MC. The continuum parameters are floated in the fit. Table U shows 
the fit results and branching fractions (BFs). The efficiencies are weighted by background-subtracted M{Kri') data 
distribution. Systematic errors are included in the likelihood by convolving the likelihood functions to calculate the 
significances and upper limits. 

5. b^d'y update by Belle 

Belle updated the measurements oi B ^ p"f, B lo^ decays using a sample of 657 million B mesons. Three signal 
modes, B^ P^l, — > and B^ oj^ are reconstructed with subdecay modes p+ — > tt+tt'^, p° tt+tt", 
uj —> TT'^TT^TT^. A helicity selection of the p and uj is applied to suppress decays with 77^/77 from B. 
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Figure 4: Projections of the fit results to AE in Afbc signal region for B+ p^, p+7 and uj^. Curves show 

the signal (dashed, red), continuum (dot-dot-dashed, blue), B — > K*"^ (dotted, magenta), other backgrounds (dash-dotted, 
green), and the total fit result (solid). 



Table II: Yield, significance with systematic uncertainty, efficiency, and branching fraction (JS) for each mode. The first and 
second errors in the yield and B are statistical and systematic, respectively. The sub-decay B{y> — > Tr^Tr^Tr") is included for 
the tJ7 mode. 



Mode 




Yield 


Significance Efficiency (%) 


B (lO"'") 


B+ ^ 
B'' ^ 


p+7 
p"7 


Q +15.2 +2.6 
U-O _i4.5 _3.9 

yc 7 +16.8 +5.1 
'J- ' -16.0 -6.1 
17 c +8.2 +1.1 
r ' -J -7.4 -1.0 


3.3 
5.0 
2.6 


8.03 ± 0.59 
14.81 ±0.95 
6.58 ± 0.76 


o 7 +2.9 +0.9 
°- ' -2.7 -1.1 
7 o +1.7 +0.9 
' -1.6 -1.0 

4.0 11:? ±1.3 



A fit to Afbc and AS (and Mi^^ for the p°7 mode) for candidates satisfying | Ai?| < 0.5 GeV and Mbc > 5.2 GeV/c^ 
is performed individually for the three ^7/0^7 signal modes and the two control modes and K'^^. Significant 

K*^^, _R'*+7 background in the ^"7 and p+7 samples are shifted in Ai? from signal peak. The shift offset of this 
background is determined from K* enriched sample from data, which is also used to determine the size using known 
kaon to pion misidentification rate. The Ai? projections of fit results are shown in Fig. [H Table HIl lists the obtained 
yields and branching fractions 0]. 

The branching fractions are combined to a single branching fraction B(B (p, w)7) = B{B^ — > p+7) = 2 x 
^B{B" p°7) = 2 X ^B{B° Lu-f), and the ratio to corresponding mode B K*"f is calculated as ^^ifj^.^ = 

0.0284 ±0.0050tES:[]i]29- This result is used to relate the ratio of CKM matrix element: \Vtd/Vts\ = 0.195t5];[!ig(exp.) ± 
0.015(th.). The isospin asymmetry Ai{B p7) = 2^B{B° p°-f)/B{B+ p+-f) - 1 is calculated to be 
Ai{B p-f) = 0.92l[j;!^^t[J;|^, which agrees with BaBar ^ (note the different definition). The direct CP-violating 
asymmetry is also measured for the first time by fitting i?+ — s- p+7 and B~ p~j events simultaneously, as 
AcpiB+ p+"f) = [iV(p"7) - 7V(p+7)]/[7V(p-7) + N{p+j)] = -0.11 ± 0.32 ± 0.09. 
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